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ABSTRACT 

The BlCEP experiment was designed specifically to search for the signature of inflationary gravitational 
waves in the polarization of the cosmic microwave background (CMB). Using a novel small-aperture refractor 
and 49 pairs of polarization-sensitive bolometers, BlCEP has completed 3 years of successful observations at 
the South Pole beginning in 2006 February. To constrain the amplitude of the inflationary B-mode polarization, 
which is expected to be at least 7 orders of magnitude fainter than the 3 K CMB intensity, precise control of 
systematic effects is essential. This paper describes the characterization of potential systematic errors for the 
BlCEP experiment, supplementing a companion paper on the initial cosmological results. Using the analysis 
pipelines for the experiment, we have simulated the impact of systematic errors on the B-mode polarization 
measurement. Guided by these simulations, we have established benchmarks for the characterization of critical 
instrumental properties including bolometer relative gains, beam mismatch, polarization orientation, telescope 
pointing, sidelobes, thermal stability, and timestream noise model. A comparison of the benchmarks with 
the measured values shows that we have characterized the instrument adequately to ensure that systematic 
errors do not limit BlCEP's 2-year results, and identifies which future refinements are likely necessary to probe 
inflationary B-mode polarization down to levels below a tensor-to-scalar ratio r = 0. 1 . 

Subject headings: cosmic background radiation — cosmology: observations — gravitational waves — infla- 
tion — instrumentation: polarimeters — telescopes 



1. INTRODUCTION 

A strong indication of an inflationary origin of the universe 
would be a detection of the curl component ("B-mode") in 
the polarizati on of the CMB arising from gravitational wave 
perturbations dDodelson et alj2 009). This primordial B-mode 
polarization is expected to peak at angular scales of —2°, 
and the magnitude of the power spectrum is described by 
the ratio r of the initial tensor-to-scalar perturbation ampli- 
tudes, a quantity directly related to the energy scale of infla- 
tion. The best published upper limit is r < 0.22 at 95% con- 
fidence, derived from WMAP CMB temperature anisotropy 
measurements at large angular scales combined with con- 
straint s from Type la supernovae and baryon acoustic oscilla- 
tions dKomatsu et al.ll2009l) . Upper limits on the B-mode po- 
larization amplitude, y/£(£+ l)C* B /27r, of —0.8 fiK rms have 
been placed by WMAP at a mul tipole moment of I— 65 and by 
OUaD at £-200, respectively dNolta et al.ll2009t iPrvke et al J 
2009). These limits are still well above the expected levels of 
confusion from either polarized Galactic foregrounds in the 
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cleanest regions of the sky or from gravitational lensing that 
converts the much brighter CMB gradient ("E-mode") polar- 
ization to B-modes at smaller angular scales. 

BlCEP (Background Imaging of Cosmic Extragalactic Po- 
larization) is an instrument designed to target the expected 
peak of the gravitational-wave signature at angular scales 
around 2°. Using proven bolometric technologies and select- 
ing the cleanest available field for observation, this instru- 
ment was designed, given sufficient observation time, to be 
capable of measuring a polarization signal of 0.08 /iK rms at 
£—100 corresponding to the BB signal expected for a tensor- 
to-scalar ratio of r = 0.1. The CMB temperature anisotropy 
has a much larger amplitude of —50 /iK rms at these angular 
scales, and imperfect rejection of it in the polarization mea- 
surement could result in a residual false signal. In addition, 
errors in polarization orientations and pointing could mix the 
— 1 /LtK Zi-mode CMB polarization signal into spurious B- 
modes. This experiment therefore requires careful instrument 
characterization and calibration to minimize systematic con- 
tamination in the polarization measurement. 

This paper supplements a companion paper that presents 
the CMB p olarization power sp ectra from the first 2 years of 
BlCEP data (Chi ang et al.l2010» . The rest of this section gives 
an overview of the instrument, the observing strategy, and the 
collected data. Section 2 describes the simulations that deter- 
mine the impact of the instrumental parameter characteriza- 
tion on the cosmological results. Section 3 describes each of 
the instrumental properties and how we calibrate them. Sec- 
tion 4 discusses how we characterize the properties of noise 
in our data and quantify the noise bias through simulations. 
We conclude by identifying what we find to be the most im- 
portant systematic uncertainties for BlCEP and discussing the 
path forward toward improving these uncertainties for future, 
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more sensitive, measurements. 

1.1. Instrument design overview 

The goal of targeting the sub-//K fi-mode polarization sig- 
nal that peaks at ~2° angular scales led to an experiment de- 
sign optimized especially for sensitivity and control of sys- 
tematic errors. The des ign of BlCEP an d th e observation 
strategy are described in Yoon et al. (2006) and Keating et al. 
(2003a). BlCEP (Figure [TJ is a compact on-axis refrac- 
tor with 49 pairs of polarization-sensitive bolometers (PSBs; 
Uones et ai1l2003l) operating in atmospheric transmission win- 
dows near the CMB peak at 100 and 150 GHz with 0.9° and 
0.6° beams, respectively (Table [TJ. We observe in two fre- 
quency bands to differentiate between the spectra of CMB 
anisotropics and sources of potential Galactic foreground con- 
tamination. The Amundsen-Scott South Pole Station, at an 
elevation of 2800 m, was chosen for its atmospheric trans- 
parency, stable weather, and constant availability of an excel- 
lent observing field on the sky. Achieving one degree resolu- 
tion at 2-3 mm wavelengths requires only a 25 cm aperture, 
which is compatible with a compact forebaffle and simple im- 
plementation of calibration measurements. 

The bolometers use neutron transmutation doped (NTD) 
germanium thermistors to measure the optical power inci- 
dent on a polarization-sensitive absorber mesh. After adjust- 
ing for relative responsivities, orthogonal PSBs within a pair 
are summed or differenced to obtain temperature or polariza- 
tion measurements. Because the orthogonal PSBs observe the 
CMB through the same optical path and atmospheric column 
with nearly identical spectral passbands, systematic contribu- 
tions to the polarization are minimized. 

TABLE 1 
Bicep Instrument Summary 

Band Center Bandwidth Beam FWHM" PSBs NET'' Per Detector 
96.0 GHz 22.3 GHz 0.93° 50 530/xKcmbv^ 

150.1GHz 39.4 GHz 0.60° 48 c 450^^1;^ 

17 Full width at half maximum, average over all the beams. 
b Noise equivalent temperature; see &\ 

1 After the first year, two of the 150 GHz pairs were converted to 220 GHz. 

The PSB layout on the focal plane (Figure© was chosen so 
that a 180° rotation about the boresight completely exchanges 
the polarization coverage on the sky. At the end of the first 
year, in 2006 November, we added prototype 220 GHz feed- 
horns in place of two of the 150 GHz ones along with the ap- 
propriate filters. We also replaced four bolometers because of 
their slow temporal response, high noise level, or poor polar- 
ization efficiency. We have omitted these and other problem- 
atic PSB pairs from CMB analysis for each observing year. 
After this refurbishment, BlCEP remained cold and operated 
without interruption until the completion of the observations 
in December 2008. 

1.2. Observing strategy 

With an instantaneous field of view spanning 18°, BlCEP 
maps an 800 deg 2 field daily by scanning the boresight in 
azimuth over a 64° range at 2.8°/s with hourly 0.25° steps 
in elevation from 55° to 60°. At each elevation step, the 
telescope completes 50 back-and-forth azimuth scans, or 100 
"half-scans," making up a "scan set." The scan speed was se- 
lected so that our target angular scales (£ ~ 30-300) appear at 



0.1-1 Hz, above a significant portion of the l/f atmospheric 
noise, while limiting motion-induced thermal fluctuations at 
the detectors. 

BlCEP operates in a 48 sidereal-hour observing cycle (Fig- 
ure©, with each cycle at one of the four fixed boresight rota- 
tion angles {-45°, 0°, 135°, 180°}. A 45° rotation about the 
boresight exchanges the polarization coverages on the sky, 
providing two independent pairings of angles, each with com- 
plete polarization coverage per overlapped sky pixel. The two 
sets rotated by 180° also help to average down systematic ef- 
fects like differential pointing. 

Each 48 -hr cycle begins with 6 hours allocated for recycling 
the refrigerator, filling liquid nitrogen (every 2 days) and liq- 
uid helium (every 4 days), and performing optical star point- 
ing calibrations along with a mount tilt measurement. The 
CMB field is completely mapped once each day in two 9-hr 
blocks, with the scan order of the upper and lower halves of 
the elevation range switched such that the azimuth ranges for 
the two days are offset. Differencing the first and second days 
of a given 48-hr observing cycle tests for potential azimuth- 
fixed contamination. Overlapping coverage of the sky from 
detectors with various polarization orientations is created by 
scanning in azimuth, stepping in elevation, and rotating the 
telescope with respect to the boresight every 2 days. 

Each scan set at a given elevation is fixed with respect to a 
given azimuth range instead of tracking the field center over 
the hour-long period. This scan strategy allows for a straight- 
forward removal of any azimuth- or scan-synchronous con- 
tamination. For each scan set and each of the two scan di- 
rections, the entire timestream is simply binned in azimuth 
to form a template signal which is subtracted from each half- 
scan. 

Relative detector responsivities are measured at the begin- 
ning and end of each 1-hr fixed-elevation scan set by fitting 
the detector response to a small chan ge in l ine-of-sight air- 
mass ("elevation nods"), described in § 13.1.21 

1.3. Collected data and observing efficiency 

The BlCEP mount was installed at the South Pole in 2005 
November, the cryostat was first cooled in December, and the 
instrument captured first astronomical light a month later. Fol- 
lowing calibration measurements and tests of the observing 
strategy, BlCEP began CMB observations in 2006 February. 
The instrument operated nearly continuously until 2008 De- 
cember, when it was decommissioned to prepare for its re- 
placement by BICEP2 on the same mount, planned for late 
2009. 

Excluding any incomplete 9-hr observing blocks, BlCEP 
acquired 180 days of CMB observations during 2006, in 
which a significant fraction of the observing season was de- 
voted to calibration measurements. The amount of CMB ob- 
servations increased to 245 days in 2007 and a similar amount 
in 2008. Although there is no evidence for Sun contamination 
during the summer, we restrict our CMB analysis to data taken 
during February-November. 

The first 2.5 months of data in 2006 were excluded from 
the current analysis because a different scan strategy was be- 
ing investigated at this time. As a coarse weather cut, we have 
excluded 9-hr blocks if the relative gains derived from ele- 
vation nods vary by more than 20% rms, averaged over the 
channels. This criterion cuts clear outliers in the distribution 
of atmospheric instability. After these cuts, 117 days in 2006 
and 226 days in 2007 remain for our baseline CMB analysis 
(Table©. 



BlCEP instrument characterization 



3 




FIG. 1. — BlCEP telescope on a 3-axis mount at its lowest elevation limit of 50°, looking out through the roof of the Dark Sector Laboratory (DSL) building 
located 800 m from the geographic South Pole. A cryostat with toroidal liquid nitrogen and liquid helium tanks encloses the entire 4 K optics, including two 
high-density polyethylene lenses and corrugated feedhorns. The polarization-sensitive bolometers are cooled with a 4 He/ 3 He/ 3 He sorption refrigerator to 250 mK. 



Furthermore, 3% of PSB pair timestreams are omitted due 
to cosmic ray hits, glitches, or >3% mismatch in the rela- 
tive gain measured at the beginning and end of each 1-hr scan 
set. Accounting for the 75% scan efficiency and the sched- 
uled calibration routines, the net CMB observing efficiency 
is 60% during the CMB observing blocks and 45% overall 
during each 2-day cycle. 



TABLE 2 
CMB Observation Summary 



Year 


PSB Pairs: used (total) 


Observing Days: 


Integration 




100 GHz 150 GHz 


used (total) 


Time" 


2006 


19 (25) 14 (24) 


117 (180) days 


4.5 xIO 6 s 


2007 


22 (25) 15 (22) 


226 (245) days 


8.8x10 s s 



° Based on 18 hours per day of CMB observation at 60% net observing efficiency. 



2. CALIBRATION GOALS AND SYSTEMATIC ERROR SIMULATION 

To process the timestreams into co-added polarization maps 
with systematic errors tolerable for our target sensitivity, an 
accurate characterization of the detectors and their beams is 
essential. Imperfections in the experiment and its character- 
ization can result in false B-mode polarization signal. Many 
of these systematic effects depend in a complex way on the 
scan strategy, so analytic estimates of the impact of an instru- 
mental uncertainty on the final power spectra serve only as 
a rough guide. Using the actual analysis pipelines, we have 



simulated the most significant instrumental uncertainties to 
establish benchmarks for how precisely each property must be 
measured. The results are summarized in TableQ] Calibration 
uncertainties that affect only the power spectrum amplitudes, 
but do not cause false polarization signals, are summarized in 
Table Hand discussed in the last paragraph of this section. 

We are most concerned with systematic errors that mix 
temperature anisotropy T and Zs-mode polarization signals 
into B-mode polarization. Instrumental properties that must 
be characterized are discussed in detail in the next section. 
The response of a PSB to radiation characterized by Stokes 
[T, Q,U], which are functions of frequency v and direction 
ft, can be modeled as: 

d(t) =K t * |n(r)+ g J dvA e F(v) J dQ P(tt) 

[T+— (Qcos2i/j + Usm2ip)]\, (1) 
1 + e J 

where i/j is the polarization orientation angle of the PSB, e 
is the cross-polarization response, is the beam function, 
F(v) is the spectral response, A e is the effective antenna area, 
g is the responsivity at Hz, n{t ) is noise, and K, is the time- 
domain impulse response associated with the detector's fre- 
quency transfer function. The temporal response function of 
each bolometer must be measured to deconvolve it from the 
raw timestream. Then the relative gain within each PSB pair 
must be determined for differencing. Since we derive rela- 
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FIG. 2. — 2007/2008 layout of the BlCEP beams with nominal locations, 
FWHM, and polarization orientations. There are six sections (separated by 
curved gray lines) having alternating "Q" or "U" PSB orientations with re- 
spect to the center. The observations are performed with the focal plane ori- 
entations of —45° , 0° , 1 35° , and 1 80° counterclockwise about the boresight, 
providing two independent and complete polarization coverages of the field 
to allow jackknife tests. (A 180° boresight rotation provides complete cover- 
age of Q and U Stokes parameters on the sky at each beam location.) 



tive gains from the atmospheric signal in elevation nods, we 
must verify that the spectral response of each PSB pair is well 
matched. In addition, the polarization differencing requires 
that the two PSBs have well matched beam shapes and point- 
ing. Finally, construction of the polarization map requires the 
knowledge of the polarization orientations of the PSBs and 
the telescope pointing. 

We established our calibration benchmarks based on Bl- 
CEP's design goal to be capable of measuring polarization 
down to levels corresponding to a tensor-to-scalar ratio r = 0.1 
without being limited by systematic effects. At r = 0. 1 , the B- 
mode polarization power spectrum, C\ b ', would have a peak 
amplitude of 1(1+ l)Cf '/2ir = 0.007 ^K 2 at I = 70-1 10. The 
benchmarks for the instrumental properties and characteriza- 
tion correspond to the values that result in spurious B-mode 
signal at the level of r = 0.1 in the simulations. 

The simulation procedure uses the same data processing 
pipelines as the main CMB power spectrum analysis, and is 
basically identical to the signal-only simulations used in de- 
termining the £-space filter function. We verified the simula- 
tion results using our two independent pipelines: one using 
QUaD's ps eudo-Q estimator on a flat sky, and the other us- 
ing Spice dChon et al.ll2004h on a curved sky. 

We begin with a ACDM model generated by 
CAMB (Lewis et al.l 120001) . using c osmological parame- 
ters derived from WMAP 5-year data dHinshawet al.ll2009t) 
and r = 0. From this model, we gene rate an ensemble o f 
simulated CMB skies using synfast dGorski et al.ll2005h . 
We then simulate observations of the BlCEP field on these 
synfast skies using pointing data from the actual scan 
patterns. We vary instrumental parameters for the PSB pairs 
with a distribution across the array corresponding to the given 



uncertainty, but constant in time. We have assumed random 
distributions, although the differential beams could have a 
pattern across the focal plane resulting in more or less false 
signal. For differential pointing, we use the actual measured 
quantities in the simulations. 

To simulate the coupling betw een non-ideal bea ms and the 
sky, we follow the formalism in iBock et all (120081) . A PSB 
timestream sample is expressed as a convolution of the beam 
with a second-order Taylor expansion of the sky signal around 
the pointing center — the first and second derivatives of simu- 
lated T, Q, and U maps are calculated with synfast. Using 
this technique, the beam convolution can be simulated quickly 
while using the exact scan trajectory of each detector, which 
is essential in quantifying the beam mismatch effects, as these 
depend on the scan strategy and on the combination of detec- 
tors with different characteristics. This formalism also allows 
simulation of the impact of uncertainties in the knowledge of 
PSB orientations and cross-polarization response. 

The simulated observation is performed at each of the 4 
boresight rotation angles and with the appropriate instrument 
configuration for each observing year. The simulated signal- 
only timestreams are fed through the pipelines to be filtered 
and co-added into maps with exactly the same weights as with 
the real data. A baseline set of spectra are computed without 
any simulated errors to (1) determine the amount of E-B leak- 
age due to timestream filtering and any other effects intrinsic 
to the pipeline, and (2) obtain the transfer functions to be ap- 
plied to the raw spectra. 

The power spectra of maps made with simulated er- 
rors are compared with those without any errors, and any 
differences — in particular, excess Z?-mode polarization power 
that has leaked from T or E — are attributed to the systematic 
errors. The simulations are performed with at least 10 realiza- 
tions of input CMB maps. 

Table [3] summarizes the instrument properties and r = 0.1 
benchmark levels for their characterization, as well as the re- 
sults of measurements described in the next section. Each 
instrumental property has been characterized to a level of pre- 
cision at least comparable to the r = 0.1 benchmark and in 
most cases much better. 

Finally, some calibration uncertainties affect only the scal- 
ing of the spectra but do not result in spurious polarization 
signals (Table|U). These effects can be calculated analytically, 
and are described in the next section. We set the benchmark 
for each of these uncertainties such that the contribution to the 
calibration of power spectrum amplitudes is AQ/Q < 10%, 
a standard easily met by our instrument characterization. 

3. INSTRUMENT CHARACTERIZATION 

We use the benchmarks found in the previous section to 
guide our program to characterize the instrumental perfor- 
mance. To reach the required level of precision, we designed 
and implemented a number of techniques for the calibration 
measurements. This section describes the measurements for 
relative gain calibrations, beam characterization, polarization 
calibration, telescope pointing, far sidelobes, and thermal sta- 
bility and compares the results with the benchmark values 
leading to a reference level of uncertainty in BB power spec- 
trum. The characterization of the instrument noise and its im- 
plications for uncertainty in the power spectrum is discussed 
in §|4] All the instrumental quantities have been measured 
with sufficient accuracy for the sensitivity achieved with the 
first 2 years of data. 
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FIG. 3. — A 48-hr observing cycle consists of 6 hours for cycling the refrigerator (not shown), two 9-hr blocks of raster scans on our main CMB field on the 
first day, a 6-hr block of similar raster scans on the Galaxy, and two more 9-hr blocks on the main CMB field on the second day. Each of the upper and lower 
halves of the elevation range is scanned in different azimuth ranges between the two days, allowing a jackknife test for ground contamination. At the beginning 
and end of each set of 50-minute azimuth scans, a ±0.6° elevation "nod" is performed to measure relative gains of every bolometer. This 48-hr cycle is repeated 
at different boresight rotation angles: {-45° ,0° , 135° , 180°}. 



TABLE 3 

Systematic errors potentially producing false B-mode polarization 



Instrument Property 


Benchmark" 


Measured 


Measurement notes 


Reference 


Relative gain uncertainty: &(g\ / gi) / (gl 1 gl) 


0.9% 


< 1.1% 


Upper limit, rms error over the array.* 


ED 


Differential pointing: ( ri — T2 )/""'' 


1.9% 


1.3% 


Average, each repeatedly characterized to 0.4% precision.'' 


EH 


Differential beam size: (<T[ -<t2)/<t 


3.6% 


< 0.3% 


Upper limit, rms over the array. 


EH 


Differential ellipticity: ( ey —e^ )/2 


1.5% 


< 0.2% 


Upper limit, rms over the array. 


EH 


Polarization orientation uncertainty: Aip 


2.3° 


< 0.7° 


Upper limit, rms absolute orientation error over the array. 


EH 


Telescope pointing uncertainty: Ab 


5' 


0.2' 


Fit residual rms in optical star pointing calibration. 


EH 


Polarized sidelobes (100, 150 GHz) 


-9, -4 dBi 


-26, -17 dBi 


Response at 30° from the beam center. 


EH 


Focal plane temperature stability: ATpp 


3nK 


1 nK 


Scan-synchronous rms fluctuation on ^~100 time scale. 


EH 


Optics temperature stability: A7pj 


4/JC 


0.7 fiK 


Scan-synchronous rms fluctuation on £~100 time scale. 


EH 



17 Benchmarks correspond to values that result in a false B-mode signal of at most r = 0.1. For /' = 0.01, all benchmarks would be lower by vlO. 

If relative gain errors are detected, we anticipate removing their effects in future analyses using a CMB temperature template map. 
c a = FWHM/^Witf2j = (0.39°, 0.26° ) at ( 100, 150) GHz. 

This measurement of differential pointing could be used in future analyses to remove the small predicted leakage of CMB temperature into polarization maps. 



TABLE 4 
Calibration uncertainties 
affecting the power spectrum amplitudes only 



Calibration Quantity 


Benchmark" 


Measured 


Absolute gain Ag/g 


5% 


2% 


Cross-polarization response Ae 


0.026 


0.01 


Relative polarization orientation A(ipi — ij)2) 


9° 


0.1° 



17 Benchmarks correspond to 10% uncertainty in the polarization power spectrum 
amplitude. 

3.1. Relative detector gains 

Polarization measurement with BlCEP relies on PSB pair 
differencing. The relative gains within each pair must be ac- 
curately determined to prevent unpolarized signal from leak- 
ing into the polarization measurement. The relative ampli- 
tudes of the CMB temperature and polarization anisotropies 
place stringent requirements on the relative gain calibration of 
the differenced detectors. Simulations of relative gain errors 
suggest that the level of false polarization depends strongly 
on how gain errors are distributed across the PSB array. With 
the best estimated distrib ution of the measured gain uncertain- 
ties described in 33.1.21 the simulations indicate that relative 



gains need to be accurate to 0.9% rms to limit the leakage 
of CMB temperature anisotropy into Z?-mode power at a level 
corresponding to r = 0. 1. 

We correct for relative gain differences in two steps, first 
by deconvolving the temporal transfer function of each PSB 
to account for frequency-dependent gains, and then by cor- 
recting for the DC gains through elevation nods. Since ele- 
vation nods use the atmospheric emission, which has a dif- 
ferent emission spectrum than the CMB temperature fluctu- 
ations, the spectral response of each PSB in a pair must be 
precisely matched. We measured the detector transfer func- 
tions each year, calibrated the DC gains hourly with elevation 
nods, and made spectral response measurements once before 
deployment and once in the field. The following subsections 
describe the measurements of temporal transfer function, DC 
responsivities, and spectral response. The transfer functions 
for each pair were measured to within 0.3% uncertainty and 
elevation nod relative gains to within 1.1% rms uncertainty 
over the PSB pairs. While the current upper limit on the pos- 
sible relative gain errors slightly exceeds the 0.9% benchmark 
for r = 0. 1 , this effect is not significant for the 2-year CMB re- 
sults; furthermore, if a significant signal was detected it could 
be corrected with a more sophisticated analysis. 
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3.1.1. Temporal transfer functions 

Analysis of the time series from each detector begins 
by deconvolving the temporal response using the measured 
frequency-domain optical transfer function of the detector. 
Since the transfer function is proportional to the gain of the 
detector as a function of frequency, it directly affects the rela- 
tive gains of a PSB pair to be differenced. The relative transfer 
functions must thus be measured with errors below the 0.9% 
benchmark set for the relative gains. 

At the nominal scan speed of 2.8°/s in azimuth at ~60° el- 
evations, our target angular scales of I = 30-300 fall into the 
frequency band of approximately 0.1-1 Hz. Since the eleva- 
tion nods described in the following section are sensitive to 
relative fluctuations at ^0.02 Hz, the transfer functions were 
measured down to 0.01 Hz. 

The primary measurement technique involved analyzing the 
step response to a fast-switched square-wave source (Gunn 
oscillator or broadband noise source) operating at 0.01 Hz, 
while under optical loading conditions representative of CMB 
observations (Figure|4]l. Possible dependence on background 
loading and detector non-linearity were explored by repeating 
the measurement with extra loading from sheets of emissive 
foam placed in the beam in combination with different signal 
strengths. The ratio of the Fourier transforms of the time- 
domain detector response and of the input square wave were 
averaged for each detector to obtain the transfer function. The 
measured transfer functions for a representative and anoma- 
lous PSBs are shown in Figure |5J they were stable against 
different combinations of loading and signal levels. 

The relative gain uncertainty due to measurement uncer- 
tainty is found to be <0.3% over the frequency range of 0.01- 
1 Hz. The measured transfer functions fit the following model 
as a function of frequency ui, 



K(uj)(x 



1-a 



(l-ioJT 1 )(l-iLUT 2 ) (l-iu>T a ) 



(2) 



where T[ 2. Q are time constants and a is the fractional amount 
of a slow additive component. However, we measured the 
transfer functions with sufficiently high signal-to-noise ratio 
so that they could be directly inverse Fourier transformed to 
define the deconvolution kernels. The median time constants 
were t\ ~ 20 ms and T2 ~ 5 ms, and r a is generally 100- 
200 ms with a typically < 0.05. From the first observing year, 
6 channels at 150 GHz were excluded from CMB analysis due 
to excessive roll off between 0.01 and 0. 1 Hz (large a and r a ). 
Two of the worst were in a single PSB pair and were replaced 
at the end of the year. 

Between the first two observing years, the transfer function 
measurements generally agreed to within 0.5% rms across the 
signal band. Two exceptions were excluded from the first year 
CMB data where the transfer functions were less well con- 
strained at lo w frequencie s. Details of the measurements and 
analysis are in Yoon (2007). 

3.1.2. Relative responsivities 

Relative gains are derived from elevation nods performed 
at the beginning and end of every 1-hr constant-elevation scan 
set (Figure|3j. The telescope scans in elevation with a rounded 
triangle wave pattern over a range of ±0.6°, varying the op- 
tical loading by ~ ±0.1 K due to the changing line-of-sight 
air mass. The bolometer responses are fit to a simple air 
mass model of atmospheric loading versus elevation, T atm oc 
csc(£X), to derive the relative gains across the array. 




FIG. 4. — Modulated source in front of the telescope aperture used for the 
measurement of transfer functions. Metal washers are embedded in transpar- 
ent polypropylene foam sheet to scatter the PIN-switched broadband noise 
source signal (from upper right) into the beam while keeping the total load- 
ing similar to that during nominal CMB observations. 
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FIG . 5 . — Measured transfer functions for a representative PSB (left) and an 
anomalous PSB (right), with error bars showing rms measurement repeata- 
bility Most of the anomalous transfer functions are well described by Equa- 
tion[2]and are repeatably measured, but are conservatively excluded from the 
initial analysis. 

The elevation nod is performed slowly over 50 s to limit 
thermal disturbances on the focal plane. During the nod, the 
diagnostic "dark" PSBs not illuminated through the feedhorns 
exhibit systematic voltage responses (also seen in the focal 
plane thermistors but not in the resistor channels) that are 
^0.4% of the typical responses of the illuminated PSBs, indi- 
cating thermal contamination at this level. To reduce the ef- 
fect of the thermally-induced signals, the two elevation nods 
for each scan set are performed in opposite patterns (up-down- 
return and down-up-return) and the average response is used. 
While the two patterns result in a small systematic difference 
in the individual gains, the PSB pair relative gains are consis- 
tent to within the measurement noise of 0.3% rms. 

PSB pair differencing is able to remove common-mode at- 
mospheric fluctuations, and the relative gains are very stable 
(Figure [6j. Even over a time scale of months, the relative 
gains are stable with ~1% rms measurement noise and ex- 
hibit no systematic variation with the optical loading. Rela- 
tive gains have also been derived from correlating timestream 
atmospheric fluctuations within PSB pairs, and although these 
have greater statistical uncertainties the results are consistent 
with the elevation nods to within ±3%. 

As an additional method to track gain variations, an in- 
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frared source supported by a foam paddle is swung into the 
be am to inject a sign al of very stable amplitude, as described 
in lYoon et ail (|2006). It produced a very repeatable (0.2% 
rms) response between the beginning and end of the 1-hr scan 
sets and also showed that the individual gains are stable with 
1% rms across the full elevation range. However, because 
of the ~3 K optical loading introduced by the swing arm 
and unknown polarization of the infrared source, the relative 
gains from the flash calibrator have not been used for pair- 
differencing. 

To quantify the level of leakage of the CMB temperature 
anisotropy into pair differences, the individual PSB pair-sum 
and pair-difference maps were cross-correlated. This analysis 
performed on the yearly maps from the first 2 years showed 
no statistically significant evidence for relative gain errors in 
the data and placed an upper limit of <1.1% rms on the an- 
gular scales of interest. There is weak evidence that a small 
subset of PSB pairs, especially at 150 GHz, have excess sum- 
difference correlations. With the full 3-year data set, we an- 
ticipate placing tighter constraints on the gain uncertainties. 
The power spectra of spurious Z?-mode polarization due to the 
best-estimate distributions of 1 . 1 % rms relative gain errors are 
plotted in Figure [TOk. Although these upper limits exceed the 
signal for r = 0.1 on some scales, especially for 150 GHz, this 
systematic effect is still well below the statistical error in the 
first 2 years of data. If measured, this leakage can be miti- 
gated by projecting out a CMB temperature anisotropy tem- 
plate from the polarization maps. 

The measured individual gains are scaled for each 1-hr 
scan set such that the mean response of the detectors in each 
frequency band is constant. Finally, the absolute gains for 
converting the detector voltage into CMB temperature units 
are derived by cross-correlating maps of CMB temperature 
anisotropy measured by BlCEP and WMAP. As described 
in detail in Chiang et al.l (120101) . the BlCEP map (in volts) 
and "BlCEP-observed" WMAP maps (in K C mb) are identi- 
cally smoothed and filtered, and we compute the absolute gain 
gi for each of the frequency bands from cross-correlations in 
multipole space. Within the multipole range of I = 56-265, 
gi is nearly flat, and the average is used as a single calibration 
factor for each frequency band. We take the standard devia- 
tion of 2% to be a conservative estimate of the absolute gain 
uncertainty in CMB temperature units. The results are consis- 
tent with those from the dielectric sheet calibrator (described 
in 33.31 ). which can provide a real time absolute calibration 
with a 10% uncertainty. 




FIG. 6. — Individual 150 GHz timestreams within a PSB pair (red and blue) 
are differenced (black) in this plot using a single relative gain fit over the 
plotted 9-hr period. For the actual CMB analysis, relative gains are updated 
for every 1 -hr scan set. 



3.1.3. Spectral response 

As described above, the relative gain calibration of a PSB 
pair is based on the relative response of the detectors to the 



change in atmospheric loading from a small nod in elevation. 
Because the temperature derivative of the CMB and the atmo- 
spheric emission have different spectral shapes, the relative 
gain chosen to match the response to elevation nods may not 
be optimal for the rejection of CMB temperature fluctuations. 

The spectral response of each channel was measured using 
two separate polarized Fourier Transform Spectrometers with 
a maximum resolution of 0.3 GHz, once in the lab and once 
in the field. Within each frequency band, the spectra were 
very similar from channel to channel (average spectra shown 
in Figure [7), and the upper limit on the expected relative gain 
errors due to spectral mismatch was roughly 1 % rms over the 
array. This current upper limit does not rule out spectral mis- 
match as a source of possible relative gain errors in the small 
subset of PSB pairs. 

In addition to the main band, we verified that there is no sig- 
nificant response at higher frequencies due to leaks in the low- 
pass filters. High-pass thick grill filters with cut-off frequen- 
cies of 165 and 255 GHz were used in front of the telescope 
aperture one at a time and the response to a chopped thermal 
source was measured. 150 GHz channels showed no sign of 
leaks beyond 255 GHz down to the noise floor at -35 dB, 
while 100 GHz channels exhibited leaks at ~ -25 dB level at 
>255 GHz. The magnitude of this small (^0.3%) leak was 
consistent between the PSBs in each pair, so the effect on rel- 
ative responsivities is expected to be negligible. 




200 



FIG. 7. — Average measured spectral response for each of BlCEP's fre- 
quency bands, normalized to unity. Overplotted are the atmospheric transmis- 
sion at the South Pole (red), the CMB spectrum (green) and its temperature 
derivative (blue). 



3.2. Beam characterization 

Mismatch in the beams of a PSB pair can result in a false 
polarization signal from unpolarized temperature fluctuations. 
Beam mismatches can also lead to the mixing between E- 
mode and B-mode polarization. The difference between two 
nearly circular beams can be decomposed into three quantities 
corresponding to monopole, dipole, and quadrupole differen- 
tials: differential beam size (ci — o^/cr, differential pointing 
(y\—Y2)/(t, and differential ellipticity {e \ —e^jl, where 
are the Gaussian beam sizes of the first and second PSBs in 
a pair, a is the average beam size, and ri 2 are the centroid 
coordinates. We define ellipticity as e = (a a — (Tt,)/(cr a + 07,), 
where a a j, are widths of the major and minor axes, respec- 
tively. Beam size and ellipticity differences are sensitive to 
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the second spatial derivative of the temperature field, while 
pointing offset is also sensitive to the temperature gradient. 

For BlCEP's focal plane layout and scan strategy, simula- 
tions show that differential beam size, pointing, and elliptic - 
ity of 3.6%, 1.9%, and 1.5% rms over the array, respectively, 
will result in spurious B-mode signal at the r = 0. 1 level. The 
effect of differential pointing was simulated using the mea- 
sured magnitude and direction of beam offsets with the ex- 
pected amount of false BB power scaling as the square of 
the magnitude. Differential beam size was simulated by in- 
troducing a random distribution of beam size differences in 
PSB pairs, while keeping the average beam size the same. 
Similarly, differential ellipticity was simulated by making the 
beam of every PSB elliptical by a small randomized amount 
such that the pairs have ellipticity differences of the given rms 
while keeping the beam sizes the same. The measurements 
described below indicate that the major axes tend to be more 
azimuthal with respect to the optical axis than radial, and the 
major axes of paired PSB beams tend to align within ^15° 
rms of each other. Therefore, the major axes in the simula- 
tions were varied around the azimuthal direction by 10° rms 
to roughly simulate the observed alignment trend within each 
pair. Analytic calculations show, and these simulations ver- 
ify, that the expected false BB power scales as the square of 
differential bea m size or ellipticity and as the fourth power of 
the beam size (Shimon et al. 2008). 

The beams were mapped by raster scanning a bright source 
at various boresight rotation angles. The far field of the tele- 
scope is about 50 m from the aperture, which permitted mea- 
surements in a high bay prior to telescope deployment as well 
as with the instrument installed at the South Pole. In the high 
bay, a thermal blackbody source was used at a 40 m distance, 
consisting of a liquid nitrogen temperature load behind chop- 
per blades covered with ambient temperature absorber. At 
the South Pole, a temporary mast was installed on the rooftop 
outside of the fixed ground screen, allowing us to position a 
source at 60° elevation 10 m away. For a truly far-field mea- 
surement, an additional mast was installed on the roof of the 
Martin A. Pomerantz Observatory (MAPO), at a distance of 
200 m, and a flat mirror was temporarily mounted on top of 
the telescope to direct the beams down to allow the observa- 
tion of the distant mast as well as low elevation astronomical 
sources (Figure [8}. The sources used included an ambient 
temperature chopper against the cold sky, a broadband noise 
source, the Moon, and Jupiter. The broadband noise source is 
an amplified thermal source that is ideal for probing low-level 
effects. 

The measured beams are well fit with a Gaussian model, 
typically resulting in 1% residuals in amplitude. The fitted 
centroids are repeatable to about 0.02°, although the accuracy 
of the absolute locations is currently limited by uncertainties 
in parallax and pointing corrections while using the flat mir- 
ror. The average measured FWHMs are 0.93° and 0.60° for 
100 and 150 GHz, respectively, about 5% smaller than pre- 
dicted from physical optics simulations. The beam widths are 
measured to ±0.5% precision and vary by ±3% across the ar- 
ray. The beams have small ellipticities of e < 1% at 100 GHz 
ande < 1.5% at 150 GHz. 

The largest beam mismatch effect is a pointing offset that 
gives rise to dipole patterns in many of the differenced beams 
(Figure [9]). The median differential pointing offset is 0.004° 
at both 100 and 150 GHz, and is on average 1.3% of the beam 
size a. The offsets were repeatable between observations of 
both the broadband noise source and the Moon to within the 




FIG. 8. — Beam mapping setup on site consisted of sources mounted on the 
top of fold-over masts. When using the mast on the MAPO building (200 m 
from the Dark Sector Laboratory), a flat mirror is mounted to direct the beams 
over the ground screen. 
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FIG. 9. — Beams for each PSB pair are normalized and differenced to pro- 
duce this composite differential beam map (in the same orientation as Fig- 
ure|2}. The overplotted lines show the fitted pointing offsets magnified by a 
factor of 100. 

measurement uncertainty of 0.4% of a. 

Simulated observations with the measured pointing off- 
sets indicate a false BB with an amplitude comparable to the 
r = 0.1 spectrum (FigureflOb). although well below the noise 
level of the initial 2-year data analysis. With the magnitude 
and direction of the pointing offsets measured precisely, the 
resulting leakage of CMB temperature gradients into polar- 
ization can be estimated and accounted for in future analysis. 
Differential beam size and ellipticity are not measured with 
significance; the measured upper limits of 0.3% and 0.2% 
rms, respectively, are negligibly small (FigureflOb.d). 

3.3. Polarization orientations and efficiencies 

To construct accurate polarization maps from PSB 
timestreams, we must know the polarization orientation an- 
gle ip and cross-polarization response e of each PSB. Ac- 
curate orientations must be used in map making to prevent 
rotation of E'-mode polarization into false B-mode. Cross- 
polarization response determines the polarization efficiency 
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(a) Relative Gain Uncertainty: <1.1% 




r=0.01 



100 GHz- 
150 GHz 

- 100x150 GHz 



(b) Differential Pointing: 1.3% 




(c) Differential Beam Size: <0.3% 



10 ' 



10 




(d) Differential Elliptioity: <0.2% 



10 - 



10 - 




(e) Polarization Orientation Uncertainty: <0.7° 



10 - 




100 
Multipole / 

FIG. 10. — Spurious BB power from simulations of measured potential 
systematic errors. Except for differential pointing, all BB estimates corre- 
spond to measured upper limits. Effects of relative gain error and differential 
pointing can be corrected for in the analysis if necessary. All the potential 
systematic un certainties are me asured to be well below the 2-year constraint 
of r < 0.72 jChiang et alSTUM . 

(1 — e)/(l+e), which affects the amplitude scaling of the power 
spectrum. We developed experimental techniques to mea- 
sure these quantities by injecting polarized radiation into the 
telescope aperture at many different angles with respect to 
the detectors. The phase and amplitude of each PSB's re- 
sponse determine ijj and e, respectively. This section discusses 
the calibration benchmarks for these quantities and describes 
three measurement techniques and their results. The absolute 
PSB orientations were measured to within ±0.7° and relative 
orientation to within ±0.1°, and e was measured to within 
±0.01. 

Angles of the PSBs can vary from their design orientations 
due to the limited mechanical tolerances with which they are 
mounted. The deviation from perfect orthogonality of a pair 
simply reduces its efficiency for polarization; however, an er- 
ror in the overall orientation of the pair can lead to rotation of 
E-modes into B-modes. With the expected fractional leakage 
being sin(2A?/>), the ~1 /iK£ -modes at 1= 100 can rotate into 
false B-modes at the r = 0. 1 level of 0.08 /iK if the orientation 
measurement is off by 2.3°. This benchmark and the expected 
scaling were verified by simulations of systematic orientation 
offset of all the PSBs. The calibration procedure was designed 
to determine the polarization orientations to within a degree. 



Another factor, though less important, is that the PSBs are 
not perfectly insensitive to polarization components orthogo- 
nal to their orientations, effectively reducing the polarization 
efficiency to (1 — e)/(l + e). To achieve 10% accuracy in the 
amplitudes of the polarization power spectra, which are pro- 
portional to (1 - e) 2 /(l + e) 2 , our goal was to measure cross- 
polarization responses e to better than ±0.026. 




FIG. 1 1. — Dielectric sheet calibrator for measuring PSB orientations con- 
sists of a beam-filling polypropylene sheet and an ambient load made of a 
highly emissive black lining, subjecting the beams to partially polarized ra- 
diation. The device is mounted on the azimuth stage, which can rotate about 
the telescope's boresight when pointed at zenith. 

The polarization orientations were measured using a rotat- 
able dielectr ic sheet (Figure ITTb . modeled after the one used 
by POLAR (lO , Delll l2002). A small partially polarized signal 
of known magnitude is created by using an 18-^m polypropy- 
lene sheet in front of the telescope aperture oriented at 45° to 
the optical axis. The sheet acts as a beam splitter transmit- 
ting most of the sky radiation but reflecting a small polarized 
fraction of the radiation from an ambient load perpendicular 
to the beam. The polarized signal is small compared to the 
unpolarized sky background so that it can provide an absolute 
responsivity calibration in optical loading conditions appro- 
priate for normal observations. The ambient load is made of a 
microwave absorber lining inside an aluminum cylinder sur- 
rounding the beam splitter. The absorber is covered with a 
1/8" thick sheet of closed cell expanded polyethylene foam 
exactly as in the forebaffie (described in ^3.5t , the combina- 
tion of which has ^95% emissivity at 100 GHz. 

We use this polarization calibrator by putting it in the place 
of the forebaffie and fixing it to the azimuth mount. With 
the telescope pointed at zenith, rotating the device with re- 
spect to the cryostat modulates the polarization signal for each 
detector while keeping the beams stationary with respect to 
the sky. The off-axis beams see complicated, but calcula- 
ble, deviations from the nominal sinusoidal modulation (Fig- 
ure ll2t . This setup produces a partial polarization of ampli- 
tude proportional to (r am b - T s ^ y ), the temperature difference 
between the ambient load and the sky loading. With an 18- 
fim polypropylene film and a typical temperatures of T. im \, = 
220 K and = 10 K, the signal amplitude is ^100 mK at 
100 GHz and ^250 mK at 150 GHz, small enough to ensure 
that the bolometer response remains linear. 

The measurements were performed several times through- 
out each observing year and produced repeatable results for 
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the individual PSB orientations with 0.1° rms. The relative 
orientation uncertainty of A(^i —ip2) = 0.1° results in negli- 
gible absolution calibration error. The PSB pairs were found 
to be orthogonal to within 0.1°, and together were within 1° 
of the design orientations shown in Figure [2] PSB orienta- 
tion measurements performed before and after the focal plane 
servicing (2006 November) show a possible discrepancy, cor- 
responding to an average of 1.0° global rotation. We therefore 
conservatively assign a <0.7° rms uncertainty in the absolute 
orientation for each year. This absolute orientation angle ac- 
curacy is sufficient for measuring r <C 0.1 (Figure [TOk). 
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FIG. 12. — Response of each PSB in a pair as a function of the dielec- 
tric sheet calibrator orientation (black), plotted over by the fits (in red/green 
dashes) with polarization orientations and responsivities as free parameters. 

The cross-polarization responses e were measured using 
two methods that also independently confirmed the absolute 
orientation measurement. One method used a rotatable wire 
grid in front of the telescope window with a chopper modu- 
lating the polarized load through a small aperture between the 
ambient absorber and the cold sky (FigurefTJll. Fitting a sinu- 
soid to the individual PSB response as a function of the wire 
grid angle gives the polarization efficiency and orientation. 
The measurements for all three years gave cross-polarization 
response values with a distribution e = 0.045 ± 0.02. One 
150 GHz bolometer was omitted from analysis for having an 
e > 0.12 and was replaced at the end of the first year. 

The other method used a modulated broadband noise source 
with a rectangular horn behind a wire grid, mounted on the 
mast 200 m away (Figurefl4l. This source was raster scanned 
by each of our beams with 18 different detector orientations 
with respect to the wire grid, fitting a 2-dimensional Gaus- 
sian to each raster. The measured cross-polarization responses 
were slightly lower with a median of e = 0.038. Based on the 
scatter against the results from the first method, we assign an 
uncertainty of Ae = 0.01, which translates to 4% uncertainty 
in the polarization power spectrum amplitudes. 

3.4. Telescope and detector pointing 

Pointing errors greater than 1% of the beam size a could 
contaminate th e B-mode spectrum at the r ~ 10~ 4 level 
JHuetalJl200l . Since the amount of spurious BB power 
scales as the square of the pointing error, the r = 0.1 bench- 
mark would correspond to 0.3er, or 5' for BlCEP. This bench- 
mark was verified to be conservative by simulating a 5' rms 
shift in boresight pointing every 10 elevation steps in one of 
our simulation pipelines and finding negligible spurious sig- 
nal compared to the r = 0.1 BB signal. An optical star point- 
ing camera is used to measure the telescope boresight point- 
ing model with uncertainties more than an order of magnitude 
smaller than required to achieve our benchmark. 




FIG. 13. — Device above the cryostat window for measuring cross- 
polarization responses and PSB orientations. The window is covered with 
a metal plate with a 2-cm Eccosorb aperture, and a 10-cm diameter wire grid 
is on a rotation stage under the circular aperture of the rectangular plate. The 
chopper modulates the load between the ambient temperature and the cold 
sky. 




FIG. 14. — Another calibration source, used on top of a mast, for measuring 
cross-polarization responses (and PSB orientations). The broadband noise 
source at 100 or 150 GHz outputs power through the rectangular feedhorn 
oriented for either vertical or horizontal polarization and through a precisely 
aligned wire grid to minimize cross-polarized signal. 



An accurate derivation of the sky coordinates from the tele- 
scope encoder readings requires a precise knowledge of the 
state of the mount, including axis tilts, encoder offsets, and 
flexure. Encoder data for the three mount axes are recorded 
synchronously with the bolometer timestreams and correc- 
tions to the raw pointing data are applied during map making 
using a pointing model. The pointing model is established 
using a compact optical star-pointing camera with a 2" reso- 
lution mounted beside the main window on top of the cryostat 
and co-aligned with the boresight rotation axis. There are 10 
dynamic parameters: AZ axis tilt magnitude and direction, 
EL axis tilt, the 3 encoder zeros, amount of telescope flexure 
oc cos(£X) and oc sin(£'L), and magnitude and direction of the 
collimation error of the pointing camera itself. A complete 
characterization of these parameters requires the observation 
of at least 20 stars. To establish a pointing model during 
the Antarctic summer, the pointing camera was designed to 
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be sensitive enough to detect magnitude +3 stars in daylight. 
For maximum contrast against the blue sky, we used a sen- 
sor with enhanced near-infrared sensitivity] and an infrared 
filter cutting off below 720 nm. We used a 100-mm diame- 
ter lens that was color-corrected and anti-reflection coated for 
720-950 nm. Its 901-mm focal length results in a small 0.5° 
field of view to keep the sky background low. Careful adjust- 
ments of the CCD camera and the mirrors reduced the optical 
camera's collimation error to 2.4'. During our first Antarc- 
tic summer season, we successfully captured 26 stars down to 
magnitude +2.9 in an elevation range of 55°-90°. 

Optical pointing calibrations were performed every two 
days during the refrigerator cycles, weather permitting, as 
well as before and after each mount re-leveling. In each run, 
the telescope is pointed at 24 stars at boresight rotation an- 
gles of -45°, 45°, and 135°, and the azimuth and elevation 
offsets required to center each star are recorded. The pointing 
data are fit to an 8-parameter model (it has not been necessary 
to fit for telescope flexure) with typical residuals of 0.2' rms. 
The pointing model has been checked by cross-correlating the 
CMB temperature anisotropy patterns between the pointing- 
corrected daily maps and the cumulative map; no systematic 
offsets or drifts are detected. 

In parallel with the optical pointing calibration, the tilt of 
the telescope mount is monitored every two days using two 
orthogonal tilt meters mounted on the azimuth stage. We ob- 
served seasonal tilt changes of up to 0.5' per month, possi- 
bly due to the building settling on the snow, and typically re- 
leveled the mount before the tilt exceeded 1'. 

Finally, to co-add maps made with different PSB pairs, the 
actual locations of all the beams relative to the boresight must 
be determined. This was accomplished by first making a full 
season co-added map of the CMB using the design locations 
and then cross-correlating the temperature anisotropy pattern 
with single detector maps for each of the four boresight rota- 
tion angles to adjust the individual beam coordinates. These 
adjusted coordinates with respect to the boresight were then 
used to iterate this process until every individual pair map was 
consistent with the full co-added map. This derivation of the 
absolute beam locations resulted in an uncertainty of 2' rms, 
based on the agreement between the first and second years. 
For BlCEP, using the CMB temperature fluctuations proved 
to be more effective than attempting a similar procedure with 
Eta Carinae, the brightest compact source accessible. 

3.5. Sidelobe rejection 

Sidelobes of the telescope beams can pick up emission from 
the bright Galactic plane and structures on the ground, possi- 
bly resulting in contamination of the polarization maps. The 
ground shields were designed to reject the ground radiation to 
a level where the contamination is below our target B-mode 
polarization sensitivity. This section describes the ground 
shield design, the sidelobe measurement, and the possible po- 
larization contamination due to the sidelobes. BlCEP's side- 
lobes are sufficiently low to enable the measurement of B- 
mode polarization to the level of r=0.01. 

BlCEP uses two levels of shielding against ground radi- 
ation: an absorptive forebaffle fixed to the cryostat and a 
large stationary reflective shield surrounding the telescope 
structure (F igure [T5") . modeled after the POLAR experiment 
dKeating et al.ll2003bh . We designed the geometry so that any 

1 Astrovid StellaCam EX with EXview HAD CCD, 
http://www.as trovid.com/prod_details. php?pid=7 




FIG. 15. — Absorptive forebaffle and the reflective ground screen. 



radiation from the ground must be diffracted at least twice be- 
fore entering the window in any telescope orientation. 

The forebaffle is an aluminum cylinder lined with a mi- 
crowave absorber to minimize reflected radiation into the tele- 
scope. It is wide enough to clear the edge pixel beams, be- 
yond which the ZotefoarrQ window (Runy arTet al J 1200*31) is 
expected to scatter <1% of the total power. At BlCEP's low- 
est nominal CMB observing elevation of 55°, the forebaffle 
is long enough to prevent radiation from sources, particularly 
the Moon, at elevations up to 27° from entering the window 
directly. The forebaffle aperture lip is rounded with a 13 cm 
radius to reduce diffraction of the diffuse beam sidelobes. 

After testing many materials for the microwave absorber, 
we chose a 10-mm thick open-cell polyurethane foam sheet 
(Eccosorb HR[_), which had the lowest measured reflectivity 
of <3% at 100 and 150 GHz when placed over a metal surface 
(150 GHz results by W. Lu & J. Ruhl 2004, private communi- 
cation). To prevent snow from accumulating in the porous 
Eccosorb foam, it is covered with a 1.6-mm thick smooth 
polyethylene foam (Volara"3), which is attached with a sili- 
cone sealant. The combined Eccosorb HR / Volara stack was 
measured to reflect ^5% of 100 GHz radiation incident at 
45°. The additional loading on the bolometers due to emis- 
sion from the forebaffle was measured to be ~1 Krj. Since 
the absorptive baffle is fixed with respect to the detectors, its 
thermal emission is expected to be stable against the modu- 
lated sky signal. 

The 2-m tall outer screen prevents the forebaffle lip from 
seeing the warm ground. The sloped aluminum surface in- 
stead reflects any diffuse sidelobes to the relatively homo- 
geneous cold sky. The 8-m top diameter is wide enough so 
that the diffracted ground radiation will never directly hit the 
window even when the telescope is at its 50° elevation lower 
limit. The edge of the outer screen is also rounded with a 
10 cm radius to reduce diffraction. 

The sidelobe response of the telescope, including the fore- 
baffle, was measured using a chopped broadband noise source 
on the mast 10m from the telescope aperture. The telescope 
was stepped in elevation up to 60° away from the source in 
0.5° increments, making one revolution about the boresight 

2 Propozote PPA30, http://zotefoams.com/pages/en/datasheets/PPA30.htm 

3 http://www.eccosorb.com/america/english/product/40/eccosorb-hr 

4 Volara; Sekisui Voltek, http://www.sekisuivoltek.com/products/volara.php 



12 



Takahashi et al. 



and back at each step to measure a radial average of the beam. 
This measurement was performed with several source attenu- 
ations down to below -50 dB to probe the far sidelobes with 
sufficient signal-to-noise ratio while also measuring the main 
beam without saturating the detector. 

Sidelobe response maps constructed from the gain-adjusted 
pair differences show that the sidelobes are up to 50% polar- 
ized (Figure [Toll. The individual sidelobe maps can be aver- 
aged over boresight rotation angle to obtain a radial profile 
(Figure ITTb. Waves polarized parallel to the absorbing fore- 
baffle lip surface are more strongly diffracted than those po- 
larized perpendicular to the lip. This difference appears to be 
mainly responsible for the polarized response in the sidelobes. 
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FIG. 16. — Map of fractional polarization for a 100-GHz feed ("CI") up 
to 80° from the beam center, showing the gain-adjusted pair difference di- 
vided by the pair sum. Beyond ^15° from the boresight where the forebaffle 
cuts off the beam, the sidelobes are generally up to 50% polarized, but with 
smooth quadrupole pattern aligned with the polarization sensitivity. 
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FIG. 17. — Azimufhally-averaged sidelobe response for the 100-GHz cen- 
tral feed. A forebaffle with an absorptive lining cuts in at 15.5° and provides 
up to an additional ~15 dB attenuation. When the telescope is at its lowest 
elevation of 50°, the lip of the outer ground screen is ~30° from the beam 
center of the central feed. 

To quantify what fraction of the total power in the beam 
remains outside of a given angle from the beam center, the 



net beam profile for horizontal and vertical polarizations was 
integrated (Figure [T8l). Less than 0.1% of the power in the 
beam is left beyond 15° and 20° of the beam center at 100 
and 150 GHz, respectively. 
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FIG. 18. — The su m of the beam profiles for the horizontal and vertical 
polarizations in Figure [T7l is integrated over the solid angle from 80° to the 
given angle, and the fraction with respect to the total power is plotted. Less 
than 0.1% of the power in the beam is left beyond 15° and 20° of the beam 
center at 100 and 150 GHz, respectively. 

Polarized sidelobes can result in spurious signals by cou- 
pling to emission from the Galaxy and the outer ground screen 
itself. To evaluate the far sidelobe rejection performance, the 
measured level of polarized response was convolved with a 
model map of a potential contaminant, and an angular power 
spectrum was computed in the BlCEP field to compare with 
the r = 0.1 BB spectrum. For the Galactic model, we used 
the dust emission predicted by Finkbeiner et al. ( 1999). The 
resulting B-mode contamination was found to be at least 400 
times below the r = 0. 1 level, meaning the measured sidelobes 
are at least 13 dB below the benchmark. The same exercise 
was repeated for a conservative model of snow accumulation 
on the ground screen panels. The contamination was even 
smaller, with the achieved rejection level at least 23 dB better 
than the benchmark. 

We have also probed potential ground contamination in our 
data through a jackknife test comparing maps made in dif- 
ferent azimuth ranges, a nd ha ve seen no evidence of ground 
signal. As mentioned in 31.21 the scan range is fixed with re- 
spect to ground during each 1-hr scan set so that subtracting 
a scan-synchronous template each hour removes any ground- 
fixed signal. 

3.6. Thermal stability 

The thermal and optical responsivities of PSBs in each pair 
are not perfectly matched, so the temperatures of the detector 
focal plane and the emissive optics must be sufficiently sta- 
ble to prevent the introduction of scan-synchronous thermal 
signals. We have measured the thermal responsivity of every 
bolometer and compared the mismatches with the focal plane 
temperature stability, which we control with a feedback loop. 
The thermal stability of both the focal plane and the optics is 
found to be adequate compared to the r = 0.1 benchmark. 

The bolometers' responsivities to the bath temperature are 
measured by correlating the detector timestreams with the 
10 mK drop when the temperature control heater is turned 
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off at the end of each refrigerator cycle. The median thermal 
responsivity, after converting voltages into CMB temperature 
units, is 0.8 /iKcmb/uKfp, and the median mismatch within 
PSB pairs is 0.08 ^Kcmb/iiKfp- Because the pair differen- 
tial responsivities are distributed randomly in the array, the 
effects of the mismatch will average out when the maps are 
co-added. The averaged mismatch over the array, consider- 
ing both the magnitudes and signs of the thermal response, is 
0.025 and 0.001 ^K C MB/nK FP at 100 and 150 GHz, respec- 
tively. To meet the r = 0.1 target of 0.08 ^K C mb at £~100, 
thermal instabilities in the focal plane must then be controlled 
to better than 3 nK rms. 

To mitigate the thermal fluctuation effects, the focal plane 
temperature is stabilized at 250 mK using a 100 kQ resis- 
tor as a control heater (nominally depositing ~0. 1 /iW) in a 
PID feedback loop with a sensitive NTD germanium thermis- 
tor. The PID parameters are set such that no active regulation 
takes place within the observational signal band of 0.1-1 Hz; 
only long time-scale drifts are controlled so that the PSB rel- 
ative gains remain unaffected. The focal plane is equipped 
with 6 pairs of monitor thermistors spaced evenly around its 
perimeter. During the first year, the thermal control scheme 
used a thermistor closest to the thermal strap connected to the 
refrigerator. To improve the recovery time from major ther- 
mal disturbances and other transient events, additional con- 
trol thermistors were installed near the control heater on the 
thermal strap at the end of the first year. The rigidity of the 
thermal straps were improved by using Vespejj supports to 
reduce susceptibility to vibrationally induced heating. Along 
with an increased response speed of the control loop, the tem- 
perature stability measured by the focal plane thermistors was 
improved from the first year to the second. 

The temperature stability of the focal plane was found to 
vary with azimuth scan speed as well as the cryostat orienta- 
tion about its axis. The stability was investigated under a va- 
riety of telescope operating conditions, including scan speeds 
in a range of 1.0°/s - 4.0°/s, and 16 evenly spaced boresight 
rotation angles. Based on the minimum variance of the scan- 
synchronous thermistor signals, we selected the 2.8°/s nomi- 
nal scan speed and four cryostat orientations about the bore- 
sight: {-45 ,0°,135°,180 }. The measured level of thermal 
fluctuations at the BlCEP focal plane was 1 nK rms in the fre- 
quency range corresponding to I = 100. 

Finally, since emission from BlCEP's optics is expected to 
be largely unpolarized, the main concern with optics temper- 
ature drifts is in mis-calibration of PSB pair optical relative 
gai ns, wh ich have an upper limit of 1.1% rms, as described 
in 33.1.21 To limit the pair-difference response to less than 
0.08 /iKcmb, the scan-synchronous optical loading fluctua- 
tions must be under ^8 /iKcmb, requiring the optics temper- 
ature to be stable to at least 4 /xKrj rms (for the 150 GHz 
band). Scan-synchronous fluctuations averaged over all the 
individual bolometers for a two month period show 0.7 /iKrj 
rms variation in the frequency range 0. 1-1 Hz. 

4. NOISE PROPERTIES AND MODELING 

Precise characterization of noise in the BlCEP data is cru- 
cial for accurately extracting the underlying CMB polariza- 
tion angular power spectrum. The detector timestreams are 
dominated by noise which must be modeled and simulated 
to subtract the noise bias from the resulting power spectra. 
Precise subtraction is critical because any misestimation re- 

5 http://www2.dupont.com/Vespel/en_US 



suits in a systematic error in the power spectrum amplitude. 
Simulating the effect of a noise misestimate on the final BB 
spectrum from the 2-year BlCEP data, we find that a ±3% 
overall error in the noise power estimate would result in a 
bias of r = ±0.1. This translates to a benchmark of 1.5% ac- 
curacy in CMB temperature units for the pair difference noise 
simulation. Also, the simulation of noise, along with that of 
signal, determines the error bars of the CMB power spectra 
and the constraints on the fi-mode polarization amplitude. To 
remove noise bias from the raw power spectra, we simulate 
signal-free, noise-only timestreams and process them with the 
same pipeline as the actual data to compute the noise power 
spectra. This section describes characterization of the noise 
properties, simulation of noise-only timestreams, and the re- 
sulting estimates of the noise bias in the final angular power 
spectrum. The noise level is consistent with expectations and 
no significant cross-correlations in noise are observed among 
pair differences, and the simulated noise has been checked for 
consistency with the actual noise level in the data. 

4.1. Properties of noise: spectra and covariance 

To simulate noise-only timestreams with the same statis- 
tical properties as the actual data, the noise properties must 
first be modeled. Because BlCEP's raw timestreams are dom- 
inated by noise, we model the noise by simply computing the 
auto and cross spectral power distributions for detector sums 
and differences. The signal-to-noise ratio in the timestreams 
is <0.2% for the PSB pair differences and 1-10% for the 
pair sums in the 0.1-1 Hz frequency range corresponding 
to I = 30-300. The significant CMB signal in the pair-sum 
timestreams is expected to introduce some error, but we per- 
mit this in the present analysis because the uncertainty in 
noise bias is expected to be much smaller than the cosmic 
variance in the temperature power spectrum. The present 
noise model accounts for correlations among all the detec- 
tors in the focal plane, but does not attempt to include cor- 
relations between half-scans or between orthogonal Fourier 
modes within each half-scan. 

PSB pair sums and differences are modeled and simulated 
after the removal of a third-order polynomial from the 20 s 
of each half-scan that are used in the CMB map making. For 
each pair sum or difference consisting of 200 points at a 10 Hz 
sampling rate, we take the Fourier transform and multiply 
by its complex conjugate for an auto power spectrum in fre- 
quency space. We multiply the Fourier transforms of different 
pair sums or differences during each half-scan, dAdg(2/ Ai/), 
to obtain complex cross spectra in units of V 2 /Hz, where Av 
= 0.05 Hz is the frequency resolution. 

For each 1-hr constant-elevation scan set, we average the 
complex spectra of 100 half-scans. By comparing the aver- 
age spectra over the first and the second halves (50 half-scans 
each) of the scan set, we checked that the noise properties 
are approximately stationary during a 1 -hr period compared 
to the uncertainties in the averaged spectra. Similarly, aver- 
age spectra for right-going and left-going scans were com- 
pared during a scan set to verify that there is no significant 
difference. To compute the noise model, we then bin the spec- 
tra into 12 logarithmically spaced frequency bands spanning 
0.05-5 Hz. The average auto spectra over all PSB pairs in 
each frequency band are plotted in Figure [19] with the de- 
tector voltages converted to CMB temperature differences. 
The NETs are consistent with expectations, and are compa- 
rable to the best achieved in other ground-based experiments 
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dRunvan et al.ll2Q03t iHinderks et al.ll2009t) . Before differenc- 
ing, the timestreams show significant atmospheric noise be- 
low 1 Hz. This is effectively rejected by pair differencing, 
although there is a hint of excess low frequency noise, espe- 
cially at 150 GHz. All auto and cross spectra for pair sums 
and pair differences are combined to form a complex noise 
covariance matrix N(/) at each of the 12 frequency bands, 
some of which are shown in Figure l20l 

For the purpose of generating the noise model, the pair sum 
and difference timestreams are gap-filled for cosmic ray hits 
and electronic glitches. This procedure is not performed in 
making the CMB maps; any detector half-scan with a glitch 
is simply excluded. Because the noise covariance matrices 
are constructed by averaging auto and cross spectra over mul- 
tiple half-scans, excluding a half-scan for a single detector 
pair sometimes causes the matrices to become non-positive- 
definite, which prevents the Cholesky decomposition neces- 
sary for the timestream simulation process described in the 
following section. Excluding a half-scan for all detectors if 
any of them contains a glitch results in data loss of up to 70% 
for the noise model calculation. We therefore fill gaps when 
possible, and reject a half-scan for all detectors if more than 
four PSB pairs display a simultaneous glitch. 
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FIG. 19. — Average noise power spectra of all PSB pairs used in the anal- 
ysis, averaged over all observing blocks during the entire first 2 years. The 
pair-sum spectra show 1// atmospheric noise, which is rejected by pair dif- 
ferencing. NETs per detector are derived from the pair-difference average 
between 0. 1 and 1 Hz. Accounting for polarization efficiencies, these corre- 
spond to an average instantaneous (i.e., single Stokes parameter) "NEQ per 
feed" of 410 and 340 fjK-Js for 100 and 150 GHz, respectively. 



4.2. Simulation of noise-only timestreams 

For each 1-hr scan set, the measured noise covariance ma- 
trix is used in the following steps to generate simulated noise 
timestreams that reflect the modeled noise correlations. 

(1) For each of the 12 frequency bands of the noise model 
spectra, we take the complex covariance matrix N(/) for sums 
and differences of all the good PSB pairs [74x74] and com- 
pute its complex Cholesky decomposition factor L(/) [lower 
triangular 74x74] such that N(/) = L(f)V(J). 

(2) For each 20-s (200-sample) half-scan and each pair 
sum or difference, we generate the positive-frequency part 
of a complex spectrum template p [74 x 100] using normally- 
distributed random numbers whose magnitude has an expec- 
tation value of unity. 

(3) For each of the 100 positive frequency bins of this vector 
of unit-spectrum templates, we multiply the Cholesky factor 
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FIG. 20. — Example correlation matrices from the BlCEP noise model. 
Each panel shows the fractional cross-correlations among PSB sums (in- 
dexes 1-37) and differences (indexes 38-74) of 37 PSB pairs used in 2007, 
with the real and imaginary components plotted in the upper-right and lower- 
left halves, respectively. The auto-spectra values along the diagonal are nor- 
malized simply by the maximum of the 74 values. The atmosphere-induced 
correlations are visible as off-diagonal structure among the pair sums. The 
imaginary component lacks substantial power, except among the pair sums 
at the lowest frequencies where the correlations due to atmosphere can be 
phase shifted depending on the relative beam locations. We define the noise 
model in 12 frequency bands spanning 0.05-5 Hz, and four of those bands 
are shown. 

from the appropriate frequency band of the noise model: y(f) 
= L(/)p(/). The resulting spectra \(f) have the same covari- 
ance as the data: 

(v(/)vt(/)) = (L(/)p(/)pt(/) L t(/)) 

= L(/)Lt(/) 

= N(/). (3) 

(4) To ensure that these v(/) approximate the spectra of real 
timestreams, the negative frequency part is set to equal the 
complex conjugate of the positive frequency part, so that the 
real part is even and the imaginary part is odd. 

(5) We take the inverse Fourier transform of \(f) to gener- 
ate 200 samples of simulated noise time series for each of the 
74 timestreams (sum and difference for the 37 PSB pairs). 

To evaluate the accuracy of the noise model and simula- 
tion, the simulated timestreams were fed back into the noise 
modeling pipeline and the resulting power spectral distribu- 
tions and covariance matrix were compared to those of the 
real data. The complex covariance was reproduced and the 
spectral amplitudes agreed to within < 1 % rms with no signif- 
icant systematic differences. 

The above procedure is used to generate 500 realizations 
of simulated noise timestreams for the entire data set. As 
with the real data, scan-synchronous templates are calculated 
and subtracted over each set of azimuth scans, and the noise 
timestreams are then co-added into maps. Example noise 
timestreams and maps are illustrated in Figure I2T1 

4.3. Noise bias in power spectra 

The noise bias (Ni) is estimated by averaging the power 
spectra from an ensemble of many simulated noise-only maps 
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FIG. 2 1 . — Real and simulated-noise timestreams for a PSB pair over four azimuth half-scans, showing that the noise qualities are accurately reproduced. The 
simulated timestreams for all the PSB pairs over the 2 years are co-added to form the noise-only maps (right). The l/f noise causes striping in the T map, while 
the Q polarization map approximates white noise. 

for the noise levels of the current 2-year data set. The actua l 
2-year constraint on r, as reported in the Chiang et al.l (120101) 
companion paper, is r = 0.03+q j'. The jackknife-derived up- 
per limit on possible misestimation of the noise bias scales 
with the noise level. Therefore, as the noise in future data re- 
leases decreases, we can expect this internal jackknife test to 
continue to allow us to place upper limits on noise misesti- 
mation (or to detect it if present) at a level corresponding to 
roughly 1/3 of the total statistical uncertainty on r, assuming 
a noise-limited BB spectrum, and less than this if a BB signal 
is detected. The final results of power spectrum estimation 
with noise bias removal and error bars based on simulations 
of signal and noise for the BlCEP 2-year data are presented in 
the CMB results paper. 

5. CONCLUSION 

BlCEP is an experiment built with a primary goal of tar- 
geting the signature of inflationary gravitational waves in the 
B-mode polarization of the CMB at angular scales near the 
expected peak around 2°. Its novel design emphasizes sim- 
plicity and systematic control, employing a carefully baffled 
compact cryogenic refractor and relying on a simple observ- 
ing strategy of azimuth-scan modulation with periodic bore- 
sight rotation. Using BlCEP's actual data analysis pipelines, 
we have identified those aspects of the experiment's instru- 
mental and noise properties which require careful control and 
characterization. We have established benchmarks for these 
quantities corresponding to the expected Z?-mode polarization 
signal for a tensor-to-scalar ratio of r = 0.1, a value several 
times smaller than the level of statistical uncertainty of the 
BlCEP 2-year result r = .03^ or r < 0.72 at 95% confi- 
dence (Chian get al.H2.Q10b . 

The instrumental characterization reported in this paper 
shows that all studied sources of potential systematic errors 
except for the pair relative gains, differential pointing, and 
possibly the noise estimation contribute to the uncertainty in 
the measurement at a level of r < 0.01. The effects which 



(Figure l22l. The noise bias in TT is 3 orders of magnitude 
smaller than the signal, and the spectrum is sample-variance 
limited. For the TE, TB, and EB cross spectra, the noise from 
each map is mostly uncorrected, so the resulting (N() are dis- 
tributed around zero. However, the noise contributes a signifi- 
cant portion of the raw EE signal and is expected to dominate 
the BB power spectrum. The error bar in the final spectrum 
is based on the scatter in the results from an ensemble of sig- 
nal+noise simulations, and for BB is dominated by noise. 

The accuracy of the noise model can be tested by compar- 
ing the spectra of the simulated noise with those of "jack- 
knife" maps, in which two maps made with complete halves 
of the data are differenced so that they are free of CMB signal 
and therefore nominally represent the noise level in the data. 
Jackknife divisions included those based on the left/right scan 
direction (shown in Figure l2"2l. azimuth range, boresight rota- 
tion angle pairs, alternating observing weeks, 2006/2007 ob- 
serving years, and focal plane detector split (more details in 
I Chiang et al. 2010). We tested for evidence of noise bias am- 
plitude misestimation using these 6 types of jackknife spectra 
for EE and BB at 100 and 150 GHz (a total of 24 spectra) by 
comparing the sum of bandpower deviations over I = 300-500 
to those from 100 realizations of noise simulations. The set of 
jackknife spectra from the actual data are found to be consis- 
tent with the simulated distributions, even in this high I range 
where the effect of noise misestimation is expected to be the 
largest due to the bias being a rapidly increasing function of I. 
Repeating this test with an intentionally introduced ±3% scal- 
ing of the noise in the simulations, we find a clearly detectable 
departure of the sum of the actual data jackknife bandpower 
deviations from those of the simulated distributions. This al- 
lows us to place an upper limit on possible misestimation of 
the noise bias at this level, at least for a uniformly scaled error 
across the I range. 

As described above, a ±3% misestimation of noise power 
would correspond to a maximum shift in our r estimate of 0. 1 
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FIG. 22. — 150 GHz power spectra of 100 realizations of simulated signal+noise and noise only, compared to the scan-direction jackknife spectra from the actual 
data. The distributions of simulated noise spectra are consistent with the data jackknife spectra, which are expected to be signal-free and a good representation 
of noise in the data. The CMB signal simulation uses the input spectra shown. Error bars in the final spectra are determined by the scatter in the signal+noise 
spectra, which is noise-dominated for BB and EB and largely cosmic variance limited for other spectra. 



were found to be controlled to this level include differential 
beam size, differential ellipticity, polarization orientation un- 
certainty, telescope pointing, sidelobes, and thermal stability. 
In addition, effects which impact the overall calibration of 
the polarization spectra, including the absolute gain, cross- 
polarization response, and relative polarization orientation, 
have been characterized well enough to ensure that calibra- 
tion uncertainty is a small fraction of our error budget. 

Of the remaining three effects, the noise is estimated with 
sufficient accuracy to limit the uncertainty contribution to 
r < 0.1, a number that will improve as more data are added. 
The differential pointing of the PSB pairs is sufficiently small 
to meet our r = 0.1 benchmark without correction, but will 
need to be taken into account to achieve lower limits on r. 
The current uncertainty in the method we use to verify our 
calibration of relative detector gains in the PSB pairs leads to 
an upper limit in possible error on r which slightly exceeds 
our r = 0. 1 benchmark. The uncertainty will improve as more 
data are included in the analysis, and a different approach may 
be needed if any significant relative gain errors are measured. 
By employing a more sophisticated analysis which allows for 



imperfectly-corrected differential gains of the PSB pairs and 
which includes the measured differential pointing, we expect 
to be able to control all studied sources of potential systematic 
errors to levels far below r = 0. 1 . 

This practical experience with BlCEP provides a guide for 
future experiments in search for the signature of inflationary 
gravitational waves in CMB polarization. 
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